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Abstract 

We calculate the parametrically dominant next-to-next-to-leading order corrections to four- 



f^ ' ferniion production e e "*"—*■ /i v^^ud -I- X at centre-of-mass energies near the W^-pair pro- 

duction threshold employing the method of unstable-particle effective theory. In total the 
correction is small, leading to a shift of 3MeV in the VF-mass measurement. We also 
discuss the implementation of realistic cuts and provide a result for the interference of 

f~^ ' single- Coulomb and soft radiative corrections that can easily be extended to include an 

OO , arbitrary number of Coulomb photons. 



1 Introduction 

Vl^-pair production at e~e~^ colliders is a key process for measuring the VF-boson mass and 
testing the non-abelian structure of the Standard Model (SM), because of its sensitivity to 
triple non-abelian gauge couplings. The total cross section has been measured at LEP2 in a 
kinematical region spanning from the T^-pair production threshold to a centre-of-mass energy of 
207 GeV [Ij with an accuracy of ~ 1% at the highest energies. The W mass has been determined 
with an error of ~ 40 MeV reconstructing the W bosons from their decay products; bounds on 
possible anomalous couplings are less stringent, and deviations from the SM predictions have 
been constrained at the per-cent level. A detailed analysis of the VF-pair production process 
will be possible at the planned International Linear Collider, where the total cross section could 
be measured at the per-mille level [2]. The precision on the W-raass determination has been 
estimated to be ~ 10 MeV by a direct reconstruction of the VF-decay products [3J, and ~ 6 
MeV by a dedicated threshold scan [Ij. Moreover, the possibility to scan the high-energy region 
would allow to measure more precisely the size of the triple non-abelian gauge couplings. 

The estimates on the PF-mass determination rely on statistics and the performance of the 
future linear collider, and they assume that the cross section can be predicted by theory to 
sufficient accuracy in order to convert its measurement into one for the W mass. For this 
reason, radiative corrections to on-shell VF-pair production have been thoroughly investigated 
in the past and are known at next-to-leading order (NLO) since the beginning of the 1980's [5j. 
However, the W bosons being unstable, a precise theoretical prediction has to be formulated 
for a final state of stable or sufficiently long-lived particles, represented by the fermion pairs 
produced by W decay, rather than for on-shell W bosons. NLO predictions for four-fermion 
production far from the VF-pair threshold region are available since some time in the double-pole 
approximation [6l[71[8] or with further simplifications [9l[T0]. Recently, a full NLO computation 
of four-fermion production has been performed in the complex-mass scheme |lltll2] without any 
kinematical approximation; moreover, a compact analytic result around the threshold region 
has been obtained in ^3j using effective field theory (EFT) methods [III [T5| [T6]. 

In particular, the analysis performed in [13] led to the following conclusions: 1) a resumma- 
tion of next-to-leading collinear logarithms from initial-state radiation is mandatory to reduce 
the error on the W mass from the threshold scan below ~ 30 MeV; 2) the NLO partonic 
cross-section calculation in the EFT approach implies a residual error of ~ 10 — 15 MeV. Al- 
though a large component of the uncertainty at point 2) can be removed using the full NLO 
four-fermion calculation |11|, I12j . the computation of the dominant higher order corrections, 
whose contribution has been estimated to be roughly ~ 5 MeV in [13], is necessary to secure 
the 6 MeV accuracy goal (3]. 

In this paper we employ EFT techniques to calculate analytically the (parametrically) 
dominant next-to-next-to-leading order (NNLO) corrections to the inclusive four-fermion pro- 
duction process e~e'^ — > ii~V^ud + X, where X stands for an arbitrary fiavour-singlet state, 
with a three- fold goal: to improve the EFT NLO calculation [13j; to derive a result which can 
be added on top of the full NLO prediction [TT], [12] ; to reduce the uncertainty on the VF-mass 
measurement below the required 5 MeV level by including a new set of higher order corrections. 

The organisation of the paper is as follows. In section [2] we outline the structure of our 
computation, review the essential features of the EFT method and identify the set of para- 
metrically leading NNLO corrections. In section [3] we describe in detail the calculation of the 
various contributions. In section [D we show the numerical impact of our result on the inclusive 



cross section and in section O we discuss the effect of realistic cuts adopting those apphed at 
LEP2 as a template. As a by-product we explain how invariant-mass cuts can be included in 
the EFT approach. Finally, section [6] contains our conclusions and two appendices collect re- 
sults related to the renormalisation of the Coulomb potential and the electromagnetic coupling, 
and the conversion from the fixed-width to the complex-mass scheme. 

2 Outline of the computation 

We consider the inclusive four-fermion production process 

e"(pi)e+(p2) -^ ii'V^uH + X, (1) 

where X denotes an arbitrary flavour-singlet state (nothing, photons, gluons, ...), in the 
kinematical regime close to the VF-pair production threshold, s = {pi + P2)^ ~ 4M^. Here 
the total cross section is dominated by the production of two resonant non-relativistic W 
bosons with virtuality of order fc^ — M^ ~ M^u^ ~ MwTy/ <C M^, where v is the non- 
relativistic velocity and M\y and T\\r are the W pole mass and decay width. We recall that 
the relation between the pole mass and the value given by experimental collaborations Mw = 
(80.403±0.029) GeV plj is given by Mw-Mw = Twl{2Mw) + 0{aU tl8|, where aew = a/ si, 
a is the fine-structure constant, s^ = sin^ 6w and 6yj denotes the weak-mixing angle. 

In the framework of unstable-particle effective field theory (EFT) |141 [T5t [T6] the total cross 
section is obtained through a re-organised loop expansion and a kinematical expansion in the 
small parameters 



s -4M^ ^2 Fv 

4M2, "'"' Mw 






All three parameters are of the same order, and in the following we will denote them as 5 
for power-counting purposes. Obviously, the EFT expansion is not equivalent to a standard 
loop expansion; in particular, being tailored to a specific kinematical region, it allows for a 
straightforward evaluation of leading higher order corrections. In order to appreciate this 
point, let us briefly review how the EFT exploits the hierarchy of scales and deals with the 
threshold region. The key observation consists in recognising four momentum scalings in the 
centre-of-mass frame, in the spirit of the method of regions |19j . 



hard : k^ ~ \k\ ~ M^y, potential : /cq ~ Mw 5, \k\ ~ M\y yd, 

(3) 
soft : A:o ~ |A?| ~ Mw S, collinear : /cq ~ Mw, k ~ M^y5, 

where k is an arbitrary loop-integration momentum. Starting from NNLO diagrams another 
mode has to be included in addition to the ones given in ([3]) that will be calleqj 'semi-soft': 

semi-soft : /cq ~ |A;| ~ Mwv^- (4) 

Hard modes are integrated out, and their contribution is encoded in appropriate matching 
coefficients. The remaining dynamical modes appear in genuine EFT loop computations and 
their different scaling properties lead to a non-standard power counting in the small expansion 
parameter 5 involving half-integer powers. 



^This is analogous to the 'soft' mode in the NRQCD hterature whose 'uhrasoft' mode corresponds to the 
'soft' mode in our conventions. 
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Figure 1: Sample diagrams in the EFT (first line) and in the full SM (second line). External 
fermionic lines are electrons and positrons; internal fermionic lines appearing in the diagrams 
in the second line represent the fi" u^ and u d doublets. 



The total cross section of ([T]) is related through unitarity to the appropriate cuts of the 
e~ e"*" forward-scattering amplitude. Figure [T]l shows a sample diagram contributing to the 
cut forward-scattering amplitude in the full Standard Model (SM). In the EFT language, 
the leading-order (LO) EFT contribution to the forward-scattering amplitude for the helicity 
configuration e7 e^, the only non- vanishing one at LO in the non-relativistic expansion, reads 
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The corresponding cut diagram is illustrated in figure [T^. Here the operators Op and O 
account at LO for the production and destruction of a pair of non-relativistic W bosons. Their 
explicit expressions are given by 



O(o) 



■K ap 



^w 



ec2,L7 " eci,L 



Vi^l Vi^l 



(6) 



where we have introduced the short-hand notation a^^V^ = a^V +a^U . Here n represents the unit 
vector for the direction of the incoming electron three-momentum pi, 7*, with i = 1,2,3, are 
the usual Dirac matrices, the subscripts on the electron fields denote the two different direction 
labels of the collinear fields and Vfj^, with i = 1,2,3, are non-relativistic spin-1 destruction 
fields for particles with electric charge ±1, whose propagators are 



feo 



|fcp 

2Mw 
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iVv 



(7) 



Loop corrections to the unstable-particle EFT scattering amplitude are then computed using 
the vertex ([6]), the propagator ^ and interactions with potential photons, related to the 
Coulomb force between the slowly-moving W bosons, (semi-) soft, and collinear photons. 

For definiteness, let us consider selected radiative corrections to the LO EFT cut diagram 
from figure [Tti, and show how the non-standard power counting arises. From the scaling 
properties ([3l), it is clear that the propagators ^ of the potential fields J7 scale as 1/6 and the 
integration measure d^ /c as 5^'^; since each vertex ([6]) contains the coupling constant aew, we 
can associate a power-counting factor a^^ 6^'^ to the LO diagram. 

We now turn to virtual Coulomb (potential) corrections, shown in figure [Da, where we 
denote the loop momentum of the Coulomb photon by kc- The corresponding SM diagram 
is depicted in figure [1^. The photon propagator scales as 1/k^ ~ 1/6 and the integration 
measure d^/cc as 5^"; including the overall factor a ~ 5 due to the photon interaction vertices 
and two additional W propagators, we obtain a power-counting factor al^6. The situation is 
quite different for the soft real-photon corrections illustrated in figure [It (see figure [Tf for the 
Standard Model counterpart). The photon propagator, parametrised by the loop momentum 
ks, scales as 1/kg ~ 1/6'^, and the integration measure d'^kg as 6'^; each collinear internal fermion 
line scales as l/{ks ■ Pi) ~ 1/6, with i = 1,2, where pi and p2 are the incoming momenta of the 
external electron and positron; the two photon vertices introduce an overall factor a ~ (5; the 
power-counting factor associated to the soft-photon diagram is a^^ 6^''^ . 

The standard loop expansion of figure [Tfe and figure [Tf treats virtual Coulomb corrections 
and real soft-photon contributions as genuine NLO effects. The re-organised EFT expansion, 
instead, shows that Coulomb-photon corrections at the VF-pair production threshold are sup- 
pressed with respect to the LO diagram by a factor 6^'"^, and can be classified as dominant 
effective N^'^LO^^"'" terms in the small expansion-parameter 6; soft-photon corrections are 
suppressed by one power of 6, and can be treated as sub-dominant NLO^^""" effects. 

The aim of this paper is to calculate corrections to the four-fermion production cross section 
which are not included in the SM NLO result, but might be relevant to the W mass analysis at 
the 5 MeV level. It is therefore natural to turn our attention to that subset of the SM NNLO 
diagrams which corresponds to EFT terms containing an extra factor 6^'"^ with respect to the 
LO result, rather than 6"^. Corrections of the same order in the EFT power counting that 
are already included in the SM LO or NLO diagrams (non-resonant contributions encoded in 
four-electron operators in the EFT; higher dimensional production operators; etc.) are not 
considered in the present paper. Those N^"L0 radiative corrections that correspond to 
SM NNLO diagrams can be readily organised in several classes: 

Mixed hard/Coulomb corrections: This class is given by diagrams with a single-Coulomb 
photon and one insertion of a hard NLO correction to the 

• production stage: Here the LO operator ^ is replaced by the NLO expression (fig- 
ure [2^). A representative diagram of the full-SM counterpart is shown in figure [2}1; 
this correction is computed in section 13. 3i 

• decay stage: A sample diagram in the standard loop expansion is shown in figure [2h- 
The implementation of this correction in the EFT discussed in section [3^71 is denoted 
by the black dot labelled (Jdecay in figure [lb. 

• propagation stage: A sample diagram in the standard loop expansion is shown in 
figure [2]F. The implementation of this correction in the EFT discussed in section 
is denoted by the black dot labelled ^residue in figure [2)3. 
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Figure 2: Sample N^'^LO diagrams in the EFT (first line) with mixed hard/Coulomb 
corrections and corresponding NNLO diagrams in the full SM (second line) . Note the insertion 
of a NLO production operator in diagram a. The same conventions of figure [T] are adopted. 



Interference of Coulomb and radiative corrections: There are two contributions in this 
class: 

• Single- Coulomb exchange and soft photons. An EFT diagram in this class is shown in 
figure [3^. A representative diagram of the full-SM counterpart is shown in figure [3li. 
This correction is computed in section 13. 2i 

• Single- Coulomb exchange and collinear photons. These corrections (see figure [3)3 
for a representative diagram in the EFT, and figure [3^ for a counterpart in the 
full SM) vanish if the electron mass is set to zero. If a finite electron mass is used 
as infrared (IR) regulator there are further contributions from soft-collinear and 
hard-collinear modes that are computed in sections 13.21 and 13. 3^ respectively. 

NLO corrections to the Coulomb potential: The relevant diagram is given by a semi-soft 
fermion bubble insertion into the Coulomb photon (see figure [3h for the EFT diagram and 
figure [3]F for the counterpart in the full SM) . This correction is computed in section 13.61 
and appendix lAl 

This concludes our survey of the dominant NNLO corrections. In the next section we will com- 
pute the different contributions identified above. At the same order in the EFT power counting 
as the above contributions we also encounter triple-Coulomb exchange that is a NNNLO cor- 
rection in the standard loop expansion. This effect can be straightforwardly included using 
EFT methods but turns out to be numerically negligible. 
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Figure 3: Sample N'^'^LO^^"'" diagrams with single-Coulomb exchange and radiative corrections 
in the EFT (first line) and corresponding NNLO diagrams in the full SM (second line). The 
same conventions of figure [1] are adopted. 



3 Evaluation of the dominant NNLO corrections 



In section [2] the dominant NNLO corrections to four-fermion production near the W-pair pro- 
duction threshold have been identified. In this section we compute all the relevant contributions. 
Subsection 13. 1 1 reviews the known result for the all-order Coulomb Green function that will enter 
subsequent calculations and extracts the triple-Coulomb exchange correction. Interference of 
single-Coulomb exchange with corrections from soft and soft-coUinear photons are computed in 
subsection 13. 2i Subsection 13.31 contains the results for interference of single- Coulomb exchange 
with hard corrections to the production stage and with hard-collinear photon exchange. These 
computations are carried out in a form that allows in principle to include all-order Coulomb 
exchange as well. This subset of NNLO corrections is combined to a finite partial result in 
subsection 13.41 and large logarithms of the electron mass are absorbed in electron structure 
functions in subsection 13.51 Radiative corrections to the single- Coulomb exchange potential 
itself are discussed in subsection 13.61 Interference effects of single-Coulomb exchange with 
corrections to the decay and propagation stages are subject of subsections 13.71 and | 



3.1 Coulomb Green function and triple-Coulomb exchange 

It is clear that Coulomb corrections will play a privileged role in our computation. In particular, 
in the following, we will derive a result which can be applied to diagrams involving an arbitrary 
number of Coulomb photons. Therefore, in order to illustrate our formalism, we review here the 
computation of pure Coulomb corrections to the four-fermion production process ([T|), given by 






Figure 4: Pure Coulomb corrections to the forward-scattering amplitude. 



the exchange of potential photons with energy ko ~ MwS and three-momentum \k\ ~ MwvS- 
Note that these corrections correspond to insertions of non-local four-boson interactions in the 
effective Lagrangian of potential non-relativistic QED [20j. They can be summed to all orders in 
perturbation theory, as shown in figureHl in terms of the zero-distance Coulomb Green function, 
the Green function Gq {fi,f2',E) of the Schrodinger operator —V^/Mw — a/r evaluated at 

r*! = r2 = 0. 

Using the representation of the Green function given in |21j and defining E = -y/i — 2M\y and 
£\Y = E + iT\Y, we can represent the Coulomb-corrected e~e^ forward-scattering amplitude 
corresponding to the diagrams in figure H] as 
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LR 






(P2 



n 



Vl|pi-)(P1 - \nh'\p2-) 4°^(0,0;£:v^), 



(8) 



where we have adopted the canonical helicity notation \p^) 



l±7^ 



u{p). The MS-subtracted 



representation for the zero-distance Coulomb Green function with a complex argument £wi 
including all-order photon exchange, reads as follows [22]: 



Gg)(0,0;fH/)^ 






£w 



Mw 



+ a 
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In 



4 Mw^w 
7~ 
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+ 7E + V' 1 



a 



2V-^vy/Mi 



w 



(9) 

Here ^e is the Euler-Mascheroni constant, // the 't Hooft unit of mass, and ^ the Euler psi- 
function. The total cross section for the left-right helicity configuration can be obtained by 
evaluating the spinor product {p2 — \rv-^^-^^\pi—){pi — \rv-^'^^^\p2—) = 16(1 — e) M^ and by 
multiplying the imaginary part of the forward-scattering amplitude ([8]) by the LO branching- 



fraction product Br 



(0) 



.(0) 



1/27. Finally, expanding the psi function appearing 



m 



in a, we get the one- and two-Coulomb photon exchange terms 
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£w 



(10) 



where the scheme-dependent real constant of ([9]) has dropped out. These three terms, shown 
in figure m are already included in the NLO EFT computation. 

Triple-Coulomb exchange arises from a NNNLO correction in the standard loop counting, 
but contributes at order N^'^LO in the EFT power-counting. The result can be obtained 
straightforwardly by expanding ([9]) up to order a^: 



Aa2l 
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Mw 
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(11) 



The correction to the hehcity-averaged cross section Acr*-"^ = A(7^^/4 directly at threshold is 
Aa'^^^{^/s = 161GeV) = O.Olfb while the effect is even smaller away from threshold. 

3.2 Soft and soft-collinear corrections 

We consider here the radiative corrections obtained including soft-photon exchange to the 
diagrams shown in figm'e HI These 0{a) contributions to the Coulomb-corrected forward- 
scattering amplitude are related to the initial-initial state interference diagrams in the effective 
theory shown in figure [SJ where the photon energy and three-momentum scale both as kQ ~ 
\k\ ~ M\y5. Analogously to the NLO calculation in [13], diagrams with a coupling of soft 
photons to Q lines cancel as can be shown using a gauge-invariance argument. The coupling of 
soft photons to collinear electrons and positrons is given by the soft-collinear effective-theory 
(SCET) Lagrangian |25 1 126 ^ 127]. and amounts to the eikonal coupling zizien^, where n^ stands 
for the direction of the four-momentum of the electron or the positron. 





Figure 5: Soft-photon corrections to the all-order Coulomb-corrected forward-scattering am- 
plitude. 



The sum of the two diagrams shown in figure can be written in close analogy with ([5]) as 

^^1f = '-^ (1 - e) 2Xs, (12) 

where Is denotes the convolution of the zero-distance Green function and a single soft-photon 
exchange correction. Including prefactors it reads 

r - o. ,.2 ~2e f d'^k G'g\0,0;£w-k^) ..,„. 

2^5 = — 327ra M^/ U / r-7—-K tt — ; tt ; r. (13) 

Here the space-time dimension is d = 4 — 2e, /i^ = ji^e^^ /{Air) and Gq is given, to all orders 
in a but in d = 4, in Q. 

The soft integral (J13p can be evaluated in two steps: first, we approximate the momenta of 
the external electron and positron introducing pi = A/iyn+ and p2 = Mwn-, where n+ and 
n_ fulfil the relations n^ = n^ = and n_|_-n_ = 2; next, we perform the A;'^ integration closing 
the integration contour in the lower complex k^ half-plane, using Cauchy theorem and Jordan 
lemma and picking up the residue of the pole at k^ = |A;| — iO. Note that all singularities of 
Gq are located in the upper k^ half plane. The result is 



^' = - U ^raT^ri) ('^■^"^)' I "' ■- ,.;. ■ (") 



The single pole in the e plane in the prefactor of (|14p is associated with the emission of 
photons collinear to the incoming electron or positron, whose mass must be neglected in the 
soft region. The finite electron mass requires the introduction of two further regions (soft- and 
hard-collinear) [13], which convert the collinear 1/e pole into a large logarithm containing the 
electron mass. The soft-coUinear correction to the forward-scattering amplitude is 






iCxSC 



i^LR = ^72^ i^ - ^) ^^SC, (15) 



where the integral Xgc, denoting the convolution of the zero-distance Green function with 
soft-collinear emission, is the same expression as (J13p . but cannot be evaluated employing 
the same approximations as for the soft corrections. For soft-collinear emission along the 
electron direction, the photon momentum satisfies n_ • k ~ Tw, ^+ • k ~ T]y ■ (m^./Mw)'^ and 
k± ~ nieTw/Mw with opposite-pointing light-like vectors n^p. In this case, we parametrise 
the external momenta according to pi = an_|_ + /3n_ and p2 ^ an- + f3n+, and fix a and 
P in terms of the kinematical variables s and ml through the relations s = 4 (a + /3)^ and 
rrig = 4a/3. The result for the soft-collinear configuration reads 



Summing up the soft and soft-collinear corrections (J12p and (jlSp and using the explicit 
results for the integrals Is and Xsc (dH) and p6]) , we obtain 



^Cx[S+SC] ^ ^cxS , aCxSC _ 327r a^^ a _ 



- ^^^^fM^\'^ 



er(l/2-e) ^ ' \m. 



x{e^^n) / dk^^ j-^^^^ . (17) 

As expected, the e pole cancels in the prefactor of (fT7|) : the infrared sensitivity of the result is 
reflected in the large logarithms ln(M^/?7ie). 

3.3 Production-vertex and hard-collinear corrections 

We turn now to the radiative corrections to the total cross section of ([T|) obtained replacing 
the LO production operator ^ with the NLO expression 
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Here the short-distance one-loop coefficients C^ ^^ and C^ ^j^ follow from the matching pro- 
cedure introduced in [15l |28]: first, the one-loop e^/^ej/j^ —>' W~W^ scattering amplitude 
is computed at LO in the non-relativistic approximation; next, the result is matched with 
the amplitude obtained with the tree-level operator in the effective theory. For our purposes, 
the one- loop coefficient C^ ^^ is irrelevant, since the e]^e^ helicity configuration does not give 



any LO contribution, and no interference between the RL and LR configurations arises. As 
explained in [13] we need to take into account only the real part of C^ [^, 
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The finite part c, 



Rec 



(l,fin) 



(l,fin) 
p,LR 



is given explicitely in Appendix B.l of [13]; numerically, we have 



'-P,LR - -10.076 for Mw = 80.377 GeV, Mz = 91.188 GeV, top-quark mass mt = 174.2 
GeV and Higgs mass Mh = 115 GeV. The hard production- vertex contributions to the 
Coulomb-corrected forward-scattering amplitude can be readily obtained appending twice the 
one-loop coefficient (fT9]l to the expression in ([8]), 



Im^gx, 






(20) 







Figure 6: Hard-collinear photon corrections to the all-order Coulomb-corrected forward- 
scattering amplitude. 

Finally, we consider the hard-collinear photon corrections illustrated in figure [H associated 



with momentum scalings k^ ~ My/ and k"^ 



~ mz 



The collinear-photon couplings arise from 



the SCET Lagrangian, and their couplings to the W bosons are encoded in collinear Wilson 
lines Wc^,^ incorporated in the production operators by replacing the electron field e^^L by 

WciBci^L and analogously for the positron field ec2,L- Hard-collinear photon corrections to the 
forward-scattering amplitude clearly have a factorised form; their contribution can be obtained 
multiplying ([8]) by the hard-collinear factor derived in [13j . 
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Summing up the production- vertex and hard-collinear photon corrections (j20p and (j2ip we 
obtain 

Im^S'""""' - I- (-^2^ +-4g'"=) = '-^^ (1 - .)ImG<?'(0.0;£„,) x 
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3.4 The total cross section 



The combination of (J17p and (j22p represents our basic expression for the higher order contri- 
butions to the ah-order Coulomb-corrected forward-scattering e^e^ amphtude in the effective 
theory. In order to better exploit the structure of the result and show how the cancellation of 
poles in the e plane takes place, we introduce a modified '+' distribution, 
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where a is an arbitrary positive real parameter. This parametrisation allows us to cast the 
integral containing the zero-distance Coulomb Green function on the right-hand side of (fH 
into 
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with f{k) = Gq {0, 0; £w — k); therefore we can prove explicitely that the sum of (fT7|) and ([22 
is free from e poles. The full correction to the total cross section ([1]) reads 
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Now, after all poles in e have cancelled, Gq {0,0]£w) can be taken to be the four-dimensional 
expression given in ([9]). The result of [13] for the NLO correction to the total cross section at 
the partonic level can be readily re-derived inserting the zero-Coulomb exchange term of ([9]) 
into (j25p . and observing that 



Im 



1 
dk-—- 

[k\a+ 



£w — k 



£w 



4£ 



W /-±i^ In -2^ -2 . (26) 

Mw \\ Mw[ \ a J \ i ^ ^ 

The new N^'^LO^^"'^ correction follows from inserting the one-Coulomb exchange term 
of ([9]) into (I25p for the total cross section; after performing explicitely the k integration. 



Im 



1 

Ra4 



In 



4:Mw{£w - k) 



we get 



. Cx[S+H] 

^"lr 



f^^ 



27 s 



Im 



Im 



In 



In^ 



£w_\ 
Mw) 



In 



Mw) 



In 



— ) 

MwJ. 



(27) 



£w_\ 
Mw) 



2 In 



£w \ ^^ f2Mw 

Mw) V "^e 






(28) 



The dependence on the a regulator has obviously cancelled; the result is finite in the e plane, 
but the infrared sensitivity of the cross section is reflected in the large logarithms ln(2Mw/'n^e)- 
In the next section we show how to absorb the large logarithms in the electron distribution 
function. 
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3.5 Initial-state radiation 

An accurate prediction for the total cross section of ([T|) requires to resum the colhnear loga- 
rithms from initial-state radiation (ISR) . Here we apply the strategy outlined in [13] , employing 
the electron structure function Fgg provided in |29t [30] . The ISR resummed cross section for 
a given helicity configuration h is 



(^h{s) 



/ dxi / dx2T^l'{xi)T^l'{x2)ah{xiX2s) 
Jo Jo 



(29) 



where the partonic cross section a/i denotes the sum of the Born result and the higher-order 
correction, ah = fjf"'^'^ + ^h^\ with logarithms of the electron mass subtracted such that 
expansion of ()29p reproduces the NLO cross section. 

Therefore, we construct Acj^^ starting from Aa^^ of (|25p and performing a 

subtraction, 



.Cx[S+H] 



AarZ'^^''^{s) 



A^r/-^\s) 



dxT];t^'\x)a2nixs) 



(30) 



where Fee is the 0{a) term in the expansion of the conventional structure function provided 
in [MIEQ], evaluated at ^/s = 2Mw, and ci^;^ is related to the all-order Coulomb-corrected 
forward-scattering amplitude ([5]) through the usual relation a^j^ = linA'2j^/{27 s). Note that 
the zero-distance Coulomb Green function appearing in ([8]) has now to be evaluated with the 
replacement E ^ E — M^(l — x). Using Tee in the limit x — > 1, 



T\l'^'\x) 



a 
2^ 



2 In 



/ 2Mw \ 



V ^e 



■) 



[1 



^«(i 



(31) 



we can write the subtraction term appearing in (j30p as 



-2 



dxrLL.(i)(:,)aC^( 



xs] 



IGvTQg^a 



2 In 



f2Mw\ 

\ rUe J 



Im|^Gg')(0,0;£:^) + 2 



dk 



gP{0,0;£w 



k) 



[k]Mw~ 



(32) 



To obtain this result, a term involving the integral Im[ Jj^ d/c G^ (0, 0; iSiy — k)/k], which is 
suppressed by a power oi £w /^^w and corresponds to hard-collinear initial-state radiation, has 
been dropped. Replacing ([32|) and (p5]) in ([30]) and choosing a = Myi/, we finally obtain 

- + ^ + Rec;^£^ j ImG^^(0,0;^H') + 2Im y^ dk 



^lr{s) 



167ra„,„a 



27sM2, 



(33) 
where the dependence on the large logarithms \n.{2M-\^ / rUe) has cancelled out. 

Using ()27p for the one-Coulomb exchange term in ([9]), we get the N^'^LO contribution 
to the total cross section to be convoluted with the electron structure functions in (1291). 



Act 



Cx[S+H] 
LR 



27s 



^ + ^ + 2RecWtlIm 



In 



£^ 



w 



■) 



MwJ 



-F2Im 



In^ 



£w_\ 
MwJ 



(34) 
The result has been cross-checked through a direct integration of the one-Coulomb exchange 
terms. 
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3.6 Corrections to the Coulomb potential 

N^' ^LO^^^ contributions also arise from semi-soft and hard corrections to the Coulomb force 
between the two Ws. The hard corrections are closely related to the renormalisation of the 
electromagnetic coupling in the Coulomb potential; they are discussed in detail in appendix El 
As a first example for semi-soft modes, consider corrections from light fermions that consist 
of the insertion of a fermion bubble into a photon propagator exchanged between two 14/^-lines 
as shown in figure[7^ but also vertex corrections (figure[7)3) and box corrections (figure[71:). For 
a potential photon the only correction at order N^' ^LO comes from the bubble insertion where 
the fermion is semi-soft, leading to the effective theory diagram shown in figure [71i. Since the 
semi-soft fermion propagator l/fiss counts as 5~^'^ and the loop measure as d'^kss ~ 5^ the 
fermion loop insertion scales as a6^ ~ 5^. Taking the additional potential photon propagator 
~ 5~^ into account this diagram is suppressed by a factor 5 compared to the pure single- 
Coulomb exchange and is therefore a N^'^LO correction to the forward-scattering amplitude. 
The triangle diagrams and box diagrams with semi-soft fermions do not contribute at order 
N^'^LO as can be seen from a power-counting argument. For instance, in the triangle diagram 
of figure [7^ the semi-soft loop momentum can be neglected compared to the VF-mass so the 
neutrino propagator contains no momentum dependence. By the same arguments as before 
the diagram is of order 5"^ compared to the single-Coulomb exchange and does not have to be 
considered. 




Tr l^eW 






I 

I 

■WV<ri^>W\M 
11/ ^N'^^ 11/ 




Figure 7: a-c: Sample diagrams in the full SM that potentially give rise to corrections from 
semi-soft fermions. Only the insertion of a semi-soft fermion bubble (a) contributes at N^' ^LO 
and leads to the EFT diagram shown in d. 



Another N^' ^LO correction potentially arises from box diagrams with exchange of two semi- 
soft photons. Example diagrams in the full theory are shown in figure [8^-c, and a representative 
diagram in the effective theory is shown in figure Ell. After interacting with a semi-soft photon 
with momentum kgs, the W propagator turns into l/{2Mwkssfi) ~ 5~^''^. Note that kssfi ~ 
M\Y\5 ^ Ty/ so that the decay width can be treated perturbatively in the semi-soft region. 
Since the photon propagators are given by 1/A;^g ~ (5~^, the semi-soft subgraph counts as 
a^ dt^kssi^ / kl^)"^! / {2kss,of' ~ a'^5~^. The potential sub-loop of the box diagram figure [Hi has 
two fi propagators so it scales as 5^''^6~'^. Therefore the contribution of these diagrams would 
again be suppressed by 5^'"^ compared to the leading order. However, at leading order in the 
non-relativistic expansion the diagrams shown in figure Eh andEb cancel each other. For this 
it is essential that the width in the propagator can be treated perturbatively in the semi- 
soft region and is not resummed in the propagator. Diagrams that contain a quartic vertex 
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like Et contribute only at higher orders because they miss at least one factor l/kgsfl from a 
VF-propagator compared to figure [8^/b. Therefore there are no contributions from semi-soft 
photons at the order we consider here. In QED this is well known from explicit calculations of 
the Coulomb potential [32|[33]. 




WW H- 






Figure 8: a-c: Sample diagrams in the full SM that potentially give rise to corrections from 
semi-soft photons. A sample diagram in the EFT is shown in d. As discussed in the text these 
diagrams do not contribute at N^'^LO. 



The only corrections to be calculated are therefore the semi-soft fermion bubble insertion 
in the single-Coulomb exchange diagram shown in figure [71i and the hard corrections discussed 
in appendix [Al For convenience we will perform the calculation in the a{Mz) scheme (see e.g. 
|3HI34]) and convert the result afterwards to the G^ scheme used in [13j. The hard corrections 
to the forward-scattering amplitude are proportional to the single Coulomb exchange, hence 



A /iNLO-Ci 



AA 



LR 



+ -^LR ^hard ' 



(35) 



where A^^^ is the unrenormalised semi-soft amplitude shown in figure[711, and ^^R ^^^ single- 
Coulomb correction to the forward-scattering amplitude included in ([8]). The hard correction 
is given in (j67p and involves the contribution of the light fermions (leptons and quarks except 
top) to the vacuum polarisation evaluated at Mz'- 



ra(Mz) 
''hard 



Ml 



f 



a 



TT 



e^^"r(e)Re 



Ml 



T\2 



r(4 - 2e) 



(36) 



Here Cf = 1 for the leptons and Cf = Nc = 3 for the quarks, and Qf is the electric charge 
of / in units of e, such that Ylf^fQ'} ~ 20/3. The total correction to the forward-scattering 
amplitude from the semi-soft bubble and the hard correction, at all orders in e, is given by 



A 4NLO-C| 
^•^LR I 



a{Mz) - 



-3e 



ECfQ} 
f 



"L«%37.. (1 _ ^) 1 r(6)r2(2-6) 



27r 



e r(4 - 2e)r(3/2 - e) 



-3e 



r(3e)r(l/2-2e)r2(l/2 + 2e) 



<"¥) 



-2e 



£w_ 



-2e 



Re 



r(4e) 






r(l/2-e)r2(l/2 + e) 



(37) 
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Expanding ([57]) in e one gets 

A ^NLO-C 

Act 



al.a' 



'LR 



\a{Mz) 



^^c,Q)[.^.{'-§^y.. 



In 



' Mw 



+ Im 



In^ 



' Mw 



(38) 

While this result is valid in the a{Mz) scheme, the numerical results in |13] were given in the 
Gfj_ scheme. According to the discussion in appendix [A] the result ()37p can be converted to the 
latter scheme by adding another term: 



NLO-C 



NLO-Ci 



,C1 



(^LR = ^(^LR \aiMz)+(>aiMz)^G^,^(^LR 



A NLO-C I 
^'^LR \a(Mz) 



2js ^c>iMz)-*G^ -tm 



In 



£^ 



w 



MwJ 



(39) 



where the single-Coulomb exchange cross-section Ao"^^ is the order a correction in (jlOp . The 

^Q is given in (j72p in appendix |Al For the same 



result for the finite conversion factor 6a(Mz)- 

input parameters as used for the hard-matching coefficient below (|19p the numerical value is 

given by 5a(Mz)-'G^ = 4.103a. 



3.7 Decay corrections 

In p4p the radiative corrections to the imaginary part of the forward-scattering amplitude are 
multiplied by the leading-order branching fraction product Br'^(M^~ -^ fj,~ 9^)31^ ' (W^ — > 
ud ) = 1/27 to extract the correction to the flavour-specific four-fermion production cross sec- 
tion. This treatment can be shown to correspond to cutting the string of fermion bubbles 
implicitly contained in the resummed propagator and selecting only those cuts that contribute 
to the desired four-fermion final state [13] . In addition to the effects described correctly by (p4|l 
there are flavour-specific radiative corrections to the decay stage contributing at N^'^LO^^"'". 
These arise from diagrams in the full theory such as shown in figure [2^, and the cuts corre- 
sponding to real emission from the final state. 

These flavour-specific corrections are included in the EFT calculation by adding the term 



Aa 



C X decay 
LR 



n(l,eui) 



.(0) 



+ 



ud 



CI 



r 



(0) 
ud 



^crYk, 



(40) 



where the one- loop electroweak corrections to the partial decay-widths, T '^p and Y £^ , can 
be found in [l3j. Note that the numerical predictions presented in [13] already include the 
2-loop QCD corrections to the hadronic decay-width multiplied by the full NLO electroweak 
cross section (including the single-Coulomb exchange). Therefore QCD corrections do not have 
to be considered in this work. 

3.8 Residue corrections to the VF-propagators 

In the proper effective-theory treatment of the computation of matching coefficients one has to 
include a factor y/2Mw{'^Rhw)~^ fo'^ each external i7 line |131 ll5j. Here w accounts for the 
normalisation of non-relativistic fields and Rhw is the hard contribution to the LSZ residue of 
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the W propagator. Introducing an expansion of the transverse self-energy around M^ and in 
the number of loops, 

Uf{k^) = M^ ^ (^" n(""'"\ (41) 



with 5 = {k — My^)/My^ and m denoting the loop order one has to one-loop accuracy [15] 

ii,-^ = i-n(i-i). (42) 

Using the on-shell scheme for field renormalisation where Rell'^'^) = 0, the residue correction 
is purely imaginary, Rj^^yy = 1 — ilmll^^'^^ = 1 -|- iVyJ /My/ with T)J the tree- level on-shell 
width. In the EFT treatment these factors reproduce the expansion of the full renormalised 
(transverse) VF-propagator 

m = .. ... ' ^^..,. = 7 '- -^ f 1 + n(^'^) + tZ7(^)) + . . . , (43) 



fc2 - M2 - YiWQ.2^ 



;r(0) 



2Mw[r,-^ + '-^ 



where the NLO contribution to the non-relativistic normalisation factor, w^"^', can be found 
in |13j , but is not needed here since in the full NLO calculation in the fixed-width or complex- 
mass scheme no non-relativistic expansion is performed. 

At order N^ "LO the residue corrections have to be included in the single- Coulomb exchange 
diagram. This corresponds to including the proper factors of \/TihW iii the NLO matching co- 
efficient Cp ^^ of the production operator in the calculation in section [331 and the WW'y vertex 
in appendix |Al In order to separate the effects corresponding to genuine loop diagrams in the 
standard loop expansion from those included in a tree-calculation using a fixed-width prescrip- 
tion, in [T3] and this paper we departed from the matching procedure sketched above and did 
not include the factors wRhw in the computation of the hard matching coefficients. Therefore 
they have to be discussed separately, paying some attention to isolate those contributions that 
are not included in a fixed-order NLO calculation in the fixed-width or complex-mass scheme. 

The factor w is solely due to the use of a non-relativistic propagator. Since no kinematic 
expansion is performed in the calculation in the fixed-width or complex-mass scheme, all cor- 
rections of this kind are already included in [IH [T2] and do not have to be considered here. 
This leaves the corrections from the derivative of the self-energy, 11^^'^-*, in (j43p . As shown in 
figure [9] these corrections contribute in two different ways to the imaginary part of the forward- 
scattering amplitude. In this figure we indicated the propagators of the non-relativistic Ws by 
the abbreviations 

The cut through the propagator in figure[9^, 2Im[l/?7] = — r^/|r/p, is interpreted as a cut 
through a self-energy insertion implicitly contained in the resummed propagator. To select the 
flavour-specific final state, the total decay width in the numerator has to be replaced by the 
(leading order) partial decay width L __ or F j [13]. This diagram therefore corresponds to 
a single-Coulomb photon exchange diagram in the full theory with the insertion of a H^ self- 
energy, expanded to first order around the mass-shell. These contributions are not included 
in the fixed-order NLO calculation in the fixed-width or complex-mass scheme where only the 
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Figure 9: Cut EFT diagrams contributing to the residue corrections to the imaginary part 
of the forward-scattering ampHtude. The O propagator with a dot indicates the insertion of 
a n*-^'^) correction factor coming from the expansion of the propagator around the pole (I43|] . 
The propagators are defined in (j44p . Analogous diagrams with insertions at the other three 
propagators are not shown. 



self-energy evaluated at Myy is resummed in the propagator [T2]. In contrast the diagram in 

figure [lb includes a cut through a propagator modified by the factor of II^^'^^ = 
The explicit expression is given by 



iv'^^/Mw 



2Im 



V 



T.(0) 

-*- w 
' Mw 



TQ 



2Mw) 



\V\ 



(45) 



This contribution corresponds to a full-theory diagram with single-Coulomb exchange and kine- 
matical correction from the expansion of a decay matrix element around the mass-shell. This 

U2r(0) 



interpretation is confirmed by inserting the one- loop result Irall^ {k'^) - 

the decay into massless fermions into the resummed propagator on the left-hand side of (j43 



k^T^^ /Mw 9 (k^) for 



Im[-iP{k)] 



-Mw 



.(0) 
w I 



k'T\];>/M, 



{k^-M^f + M^ 



w 



(fc2r[J)/M2, 



+ 0(1). 



(46) 



A fixed-width prescription as employed in the complex-mass scheme corresponds to replacing 
k'^TyJ /M"^ by T\J in the denominator, but not in the numerator, where the factor of k"^ arises 
from the integration over the two-particle phase space of the W decay products. Decomposing 
the VF-momentum as A; = Mwv + r with v = (1,0) and a potential residual momentum r 
shows that the 0{5) correction from expanding the factor /c^ in the numerator is indeed given 
by (j45p (up to the usual normalisation factor 2Mw)- Since in the fixed-order calculation of 
four-fermion production the decays are treated without kinematical approximations, this term 
does not correspond to a genuine NNLO contribution and must not be included here. 

The only correction that has to be included therefore are four cut diagrams of the form of 
figure [9^ with the residue correction inserted at the different PF-propagators. In the following 
we use again the full on-shell width Tw rather than the tree-expression '^w- Shifting the loop 
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momenta in some of the diagrams, the sum of the four terms can be brought to the form 
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(47) 



.(0) 



Here GqI denotes the single-Coulomb exchange term in the Coulomb Green function ([9]). Terms 
that vanish upon performing the go and rg integrations by closing the integration contour in 
the upper half-plane have been dropped. The explicit computation of (l47ll leads to 



Aa 



Cxres 
LR 



27s ii% " 



2\£w\{Re£w + \£w\) 



r2 



(48) 



4 Numerical analysis 

We can now assemble the N^'^LO contribution to the total cross section of the scattering process 
e~e~^ -^ ^~u^ud + X, given by the sum of the several corrections computed in section [3], 



a 



(3/2) 
LR 



A -Cx[S+H] A NLO-C , A ^Cx decay , a ^Cxres , A ^C3 
^^L/J +^'^LR +^(^LR +^^LR +^'^LR: 



(49) 



which is to be inserted into the convolution with the electron structure functions in ()29p . Recall 
that this refers to the e^e^ helicity state while there are no radiative corrections to the other 
helicity combinations from NNLO SM diagrams that contribute at N'^'^LO in the EFT power 
counting. As input parameters for the masses of the gauge bosons and the W width we use 



Mw = 80.377 GeV, 



w 



2.09201 GeV, 



Mz = 91. 188 GeV. 



(50) 



In addition, we employ for the masses of the electron, the top quark and the Higgs boson 



■nif.. 



0.51099892 MeV, 



mt 



174.2 GeV, 



M// = 115GeV, 



(51) 



and we extract the fine-structure constant a from the relation a = \/2G^-/Vf^s^/7r, where the 
Fermi-coupling constant is G^ = 1.16637 • 10~^ GeV~ and the cosine of the weak-mixing angle 
reads c^ = Myz/Mz- 

In table [1] we show the numerical results for the individual N^'^LO corrections to the unpo- 
larised cross section, Act* = Ao"}^^/4, at various centre-of-mass energies near the Vl^-pair pro- 
duction threshold. In detail, the table contains corrections from interference of single-Coulomb 
exchange and soft /hard corrections from (j34l) (third column), from the NLO corrections to the 
Coulomb-potential (j39p (fourth column), from interference of single- Coulomb exchange with 
decay and residue correction in equations (j40p and (|48p (fifth and sixth columns) and from the 
triple-Coulomb correction in (jlip (seventh column). The result (I49p for the sum of the different 
contributions is given in the second column. It is also instructive to compare to the size of the 
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a{e e+ ^ ^ u^ud X){Po) 


Vi [GeV] 


^(3/2) 


^^Cx[S+H] 


A^NLO-c 


Ao-Cx decay 


^^Cxrcs 


AfjCs 


AaC2 


158 


-0.001 


-0.116 


0.104 


-0.037 


0.044 


0.004 


0.151 


161 


0.147 


-0.321 


0.226 


-0.091 


0.324 


0.010 


0.437 


164 


0.811 


-0.417 


0.393 


-0.134 


0.965 


0.003 


0.399 


167 


1.287 


-0.389 


0.473 


-0.142 


1.345 


0.001 


0.303 


170 


1.577 


-0.354 


0.511 


-0.142 


1.561 


0.000 


0.246 



Table 1: Combined N'^^LO corrections (second column) and separate contributions from in- 
terference of single-Coulomb exchange with soft and hard corrections (third colum), renor- 
malisation of the Coulomb potential (fourth column), interference of decay correction and 
single-Coulomb exchange (fifth column), interference of residue correction and single-Coulomb 
exchange (sixth column) and triple-Coulomb exchange (all corrections are without ISR im- 
provement). For comparison the NLO contribution from double-Coulomb exchange (C2, second 
column) are also shown. 





a{e e^ ^ ^ v^udX){'ib) 


V^ [GeV] 


Born 


Born (ISR) 


NLO 


^(3/2) 


^(3/2) 


158 


61.67(2) 


45.64(2) 

[-26.0%] 


49.19(2) 

[-20.2%] 


-0.001 

[-0.0%o] 


0.000 

[+0.0%o] 


161 


154.19(6) 


108.60(4) 

[-29.6%] 


117.81(5) 

[-23.6%] 


0.147 

[+1.0%o] 


0.087 

[+0.6%o] 


164 


303.0(1) 


219.7(1) 

[-27.5%] 


234.9(1) 

[-22.5%] 


0.811 

[+2.7%o] 


0.544 

[+1.8%o] 


167 


408.8(2) 


310.2(1) 

[-24.1%] 


328.2(1) 
[-19.7%] 


1.287 
[+3.1%o] 


0.936 

[+2.3%o] 


170 


481.7(2) 


378.4(2) 
[-21.4%] 


398.0(2) 

[-17.4%] 


1.577 
[+3.3%o] 


1.207 

[+2.5%o] 



Table 2: Two implementations of the N^'^LO corrections, which differ by the treatment of 
initial-state radiation compared to the "exact" Born cross section without (second column) 
and with (third column) ISR improvement and the NLO EFT result including ISR (fourth 
column). The relative correction in brackets is given with respect to the Born cross section in 
the second column. 
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Figure 10: Corrections relative to the Born cross section: Combined N^'^LO (solid/black), 
interference of Coulomb with soft and hard corrections (long-dashed/blue), correction to 
the Coulomb potential (dash-dotted/red), interference of Coulomb and decay corrections 
(short-dashed/green) and interference of residue correction and single-Coulomb exchange (dot- 
ted/magenta). For comparison the NLO correction from double-Coulomb exchange (light 
solid/cyan) is also shown. 



pure two-Coulomb exchange corrections given by the last term in (|10p (eighth column) that 
also arise from NNLO diagrams in the full theory but appear already at NLO in the EFT power 
counting. Figure [10] shows the size of the individual corrections and the combined corrections 
relative to the Born cross section in the full theory computed using a fixed-width prescription. 
It can be seen that individual corrections are comparable in magnitude to the second Coulomb 
correction but cancel to a certain extent, in particular in the immediate threshold region where 
the total N^' ^LO correction is about one per-mille of the Born cross section. Above threshold 
the residue correction dominates and the total correction rises to about three per-mille of the 
Born cross section. 

The final results for the N^'^LO corrections including ISR improvement (i.e. the result of 
inserting a^^'"^' into the convolution with the electron structure functions in ()29p ) is given in 
the last column of table O together with the previous result without ISR improvement shown 
in the fifth column. For comparison the Born cross section in the fixed-width prescription 
with and without ISR improvement generated using WHIZARD [35] and the result of the 
NLO calculation of [13] (including ISR improvement) are also shown. The ISR improvement is 
seen to reduce the N^' ^LO corrections by about 40% at threshold and 25% at 170 GeV. Note 
that the effect is bigger than in the case of the NLO corrections [13j that are reduced by ISR 
improvement by about 20% near threshold while there is almost no effect at 170 GeV. We 
remind the reader that when the N^'^^LO terms presented here are added to the full 1-loop 
calculation in the fixed-width or complex-mass scheme j ll^ [T2] (as shown in appendix [Bl the 
N^"L0 result can be added to the fixed-order NLO result in the complex-mass scheme without 
modification), the double-Coulomb exchange terms C2 shown in table [2] must also be added, 
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since these two- loop virtual effects, although part of the NLO EFT calculation shown in table [21 
are not included in the fixed-order NLO calculation |11|,I12J. 

In [13] the impact of the interference of single-Coulomb exchange with hard and soft cor- 
rections on the W-mass measurement was estimated as [iJMiy] ~ —5 MeV. This expectation 
was based on a naive estimate for these contributions in eq. (91) of [13] that corresponds to 
a correction to the cross section of Aug^ (IGlGeV) ~ — 0.27fb. From table [1] one sees that 
this correctly captures the order of magnitude of the contribution from interference with soft 
corrections and hard corrections to the production operator (third column), that is, however, 
almost completely cancelled by the correction to the Coulomb potential. Adopting the same 
procedure as in section 6.4 of [13] to estimate the shift of the W mass, however assigning a 
relative error to each energy point that scales as one over the square root of the expected num- 
ber of events, we find that the impact of the N^^LO corrections on the VF-mass measurement 
is about 3 MeV (5 MeV if the N^'^LO correction is not convoluted with ISR), smaller than 
the targeted accuracy. Since other SM NNLO terms are expected to be even smaller, we may 
conclude that the (partonic) four-fermion cross section near the W^-pair production threshold 
is known with sufficient precision. 

5 Towards the experimentally measured cross section 

In measurements of the four-fermion production cross section, cuts on momenta or angles of the 
observed particles have to be applied in order to disentangle signal and background processes. 
While presently we do not aim at a general treatment of such phase-space cuts using effective- 
theory methods, we would like to estimate their effect on the corrections calculated in this 
article. 

For orientation, we consider the selection cuts used for the measurement of the four-fermion 
production cross section at ^/s = 161 GeV at LEP [361 EZl EHl EH]. To be definite, we use the 
cuts used by the L3 collaboration in [38] to select the qqni'{'^) final state that can be summarised 
as follows: 

(i) The muon momentum has to satisfy |p^| > 20 GeV; 

(ii) the jet-jet invariant mass Mjj and the invariant mass M^y of the muon-neutrino system 
have to satisfy 40 GeV < Mjj < 120 GeV and M^,^ > 55 GeV, respectively; 

(iii) the angle between muons and both hadronic jets must satisfy 6^j > 15 degrees to suppress 
backgrounds from qq{'y) production where muons arise as decay products of hadrons; 

(iv) the polar angle of the missing-momentum vector has to satisfy | cos 9^1 < 0.95 to suppress 
qqil) events where the missing energy arises from a photon lost in the beam pipe. 

To discuss the possible effects of these cuts on the N^/^LO corrections, it has to be understood 
to which extent they can be incorporated in the effective field theory framework. This will be 
discussed first for the leading-order cross section before the effect on the radiative corrections 
is considered. 

We have studied the effect of the cuts on the Born cross section in the full SM using 
WHIZARD [35]. The numerical effects are shown in table [3l These results do not include ISR 
improvement but we have checked that adding ISR does not change the relative impact of the 
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Cut 


0-Born(e 6+ ^ /U U^ud){Po) 


O'cut/o'tot 


- 


154.18(5) 




Ip^I > 20 GeV 


153.71(5) 


99.69(5) % 


M^^ > 55 GeV, 40 GeV < Mjj < 120 GeV 


150.61(5) 


97.68(5) % 


d^j > 15 degrees 


149.35(5) 


96.87(5) % 


1 cos 61^1 <0.95 


148.28(5) 


96.17(5) % 


all 


140.03(5) 


90.82(5) % 



Table 3: Effects of the phase-space cuts used in [38j on the Born cross section in the full SM 
at ^/s = 161 GeV without ISR improvement computed using WHIZARD. 



cuts. Although individual cuts are not very restrictive, their combined effect reduces the cross 
section by about 9 percent. The effect of the cut ^ is of the order of three per-mille and will not 
be considered further. The remaining cuts fall into two categories: the cut (jn]) on the invariant 
masses of pairs of decay products can be implemented in the effective-theory calculation as 
discussed in subsection 15. li As shown there, these cuts do not affect the corrections calculated 
in this article at all. The cuts (jiii|) and (jiv|) are sensitive to angular distributions of the decay 
products of the W bosons and are more problematic in the approach followed in [13| and 
here. The uncertainty on the theoretical prediction introduced by these cuts is estimated in 
subsection 15.21 



5.1 Invariant-mass cuts in the effective theory 

The precise treatment of invariant-mass cuts on the VF-decay products of the form — A^ < 
M? r, — M^, < A| in the effective theory depends on the scaling assigned to the ratio A/M\y 
with respect to the expansion parameter 5 ~ Tyy/Myy. For the very loose cuts in (jn]) above, it 
is appropriate to count A/Mw ~ 1 but we find it illuminating to consider also the possibility 
of tighter cuts of the order A/Mw ~ '\/Tw/Mw ~ VsE 

To discuss the correct treatment of invariant mass cuts, recall that the total cross section is 
extracted from appropriate unitarity cuts of the e~e'^ forward-scattering amplitude [13j. The 
relevant diagrams receive contributions from VF-bosons with potential and hard loop momenta. 
The potential region contributes starting at LO and is reproduced in the effective theory by the 
matrix-element of production/destruction operators as shown for leading order in ([5]), where 
the VF-propagators are given by ([7]). As long as no soft-photon corrections are considered, the 
cut on the invariant mass of the decay products translates to cuts on the momenta of the Ws 
circulating in the loop of the form — A^ < p^r — My^ < A^. These are implemented by inserting 
a product of step- functions ^(A^ — p^r + M^) 0(Af + p^ — M^) in the cut loop integral. 
We briefly comment on soft photons below. The hard region contributes from N^' ^LO and is 
reproduced in the effective theory by four-electron operators [13] . The effect of invariant mass- 
cuts in the hard region is included by taking the step functions into account in the calculation 
of the matching coefficients of the four-electron operators. We now discuss the two different 
possible scalings assigned to A/M\y in turn. 



^ In this respect it is interesting to mention that in the effective-theory treatment of top-pair production near 
threshold [40] the assumed hierarchy F ^ A <C M is closer to the second case. 
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Loose cuts: A ~ M]y Consider the contribution of the potential region first. Decomposing 
the T^-loop momentum as pw = My^v + r with u = (1, ) and a residual potential momentum 
r = (r^ ^r) ~ [5, yd), the step functions implementing the cuts are expanded as 0{h? it 
(2Mvi/?'o — ^))- Since by assumption A » ro,r^ ~ 5 the momentum can be dropped in 
the step function so that at leading order the cut can be neglected in the loop integrals in the 
effective theory. In the calculation of the matching coefficients of the four-electron operators, 
on the other hand, the step functions are operative since the VF-loop momenta are hard and 
taken to be of the order p^ ~ M^ ~ A^. Therefore a loose cut with A ~ My/ is taken into 
account entirely by modifying the matching coefficient of the four-electron operator, while loop 
integrals in the effective theory are to be performed without constraint on the 1^-momenta. 

Tight cuts: A ~ My/y5 Here the situation is reversed compared to the previous case: in the 
evaluation of cut loop integrals in the effective theory one has by assumption \2M\YrQ — f^| ~ 
A^, so the step functions have to be taken into account. In the calculation of the matching 
coefficients of the four-electron operators, the W momenta are taken to be hard and satisfy 
Ip^-M^I > A^. Since6i(A2-|p2^-M^|) ~ 6i(-|p^-M^|) = 0, tight cuts lead to vanishing 
matching coefficients. Therefore the invariant-mass cuts have to be taken into account in 
the loop calculations in the effective theory, while in the case of tight cuts the four-electron 
operators do not contribute to the cross section at all. Since A^ ~ MyrVy/ some recombination 
prescription of soft photons and final state fermions has to be supplied. 

To verify that the above procedure is the appropriate treatment of invariant-mass cuts in 
the effective theory, we have computed the full tree-level cross section for the process ([1]) with 
symmetric invariant-mass cuts |M^j- M^| < A^, |M^p^ - M^| < A^ using WHIZARD, and 
compared it to the effective theory cross section obtained using either the counting A ~ Myr 
or A ~ My/y/5. The result shown in figure [TT] shows good agreement with the effective-theory 
treatment in the regions where the respective counting rule is appropriate. Implementing 
the cuts of [38] quoted in point ([uj) at the beginning of the section by modifying the matching 
coefficient of the four-electron operator we obtain crEFT(161GeV) = 150.71 fb in good agreement 
with the WHIZARD result of o-(161GeV) = 150.61 ± 0.05 fb. 

We can now address the effect of cuts of the form (ju]) on radiative corrections. The discussion 
above shows that at the order considered here the soft-photon and Coulomb-photon corrections 
computed in section [3] do not have to be modified by invariant mass cuts as long as they are of 
the order A ~ My/. The effects of such cuts on radiative corrections ffist appears in N^'^LO 
contributions to the four-fermion operators which we have not considered in this article since 
they are included in the NLO four-fermion calculation of |1H [T2] . It is nevertheless instructive 
to consider the expected effect using the estimate of the correction to the total cross section [13] 



^(^fPiVs = leiGeV) ~ a'^/(108s^M^) = 0.74 fb, which is consistent with the difference of 
the EFT NLO calculation and \\.\\ [T2] . Estimating the effect of the invariant-mass cuts on this 
contribution to be of the same order of magnitude of the one on the leading-order four-fermion 
operator, we have 

Aaf/' ) 1 ^ 1 ~ 0.03 fb . (52) 

Clearly the effect of the cut is totally negligible. 
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Figure 11: Comparison of the Born cross section in the full SM computed with WHIZARD 
(red dots) with the effective-theory result for loose-cut implementation (dashed blue curve) 
and tight-cut implementation (solid black curve). 



5.2 Estimating the effect of angular cuts 

To discuss the treatment of cuts that involve angles between final state particles such as ([ui|) 
and (|iv|) in the effective-theory framework we will first consider the effect of a cut on the 
angle between a final-state particle and the beam pipe such as (|iv|) on the leading-order EFT 
description before we estimate the effect on the corrections computed in the present paper. 

To discuss the effect of angular cuts we can use a variant of the effective theory that includes 
the decay products [28]. In the leading-order effective-theory description the W bosons are 
produced in a s-wave by the operator ([6]). At this order the angular distribution of the Ws 
is predicted to be isotropic. Due to spin correlations the decay-product distributions are not 
isotropic, and an explicit calculation gives (we consider for definiteness the neutrino angular 
distribution) 

^ '^"""' " '^ cos0^)(3 + cose^). (53) 



^(0) dcosi 

"EFT 



16^ 



Computing the effect of the cut ([iv]) using ([53]) gives a result k = a 



(0) 



t/^: 



(0) 



96.16% which 

is in very good agreement with the number given in table [31 As another example, consider the 
neutrino forward-backward asymmetry \{a{cos6u > 0) — a{cos9i, < 0)|/<Ttot = 3/8 computed 
from (f53]l as compared to the result ~ 30% obtained using WHIZARD. 

A closer look at the angular distributions generated by WHIZARD reveals some deviations 
from the theoretical leading-order results, for example a forward-backward asymmetry of the 
T^-bosons \{a{cos6w > 0) — a{cos6w < O)|/o"tot ~ 30%. These deviations are mainly to 
be attributed to non-resonant momentum configurations that are incorporated in the four- 
fermion operators that are formally of order N^' ^LO. At threshold the effect of the four-fermion 
operators is about 40 percent of the leading order EFT contributions [T3], which is consistent 
with the order of magnitude of the observed asymmetries. The WHIZARD results show an 
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isotropic distribution of the Ws if at the same time invariant mass-cuts in the 'tight-cut' regime 
are apphed where the four-electron operators vanish, as shown in section 15.11 

Turning to the corrections computed in this paper, we first note that the Coulomb correction 
neither modifies the angular distribution nor the polarization of the W bosons. The same 
is true for the one-loop corrections to the production operator C^ [^ and the hard-collinear 



corrections considered in section 13.31 that can be both factorized from the Coulomb Green 
function. Therefore for all corrections considered in this paper apart from the soft corrections 
computed in section 13.21 the effect of the angular cuts will be the same as for the leading-order 
result. Since the numerical effect of the N^' ^LO corrections on the total cross section is in the 
one to three per-mille range, the effect of the angular cuts (which reduce the cross section by 
about 7 percent) on these corrections will be completely negligible. 

We expect that this conclusion is not significantly changed by the soft corrections. First, 
the components of the momentum of the soft emitted photon are 0{Tw) (see ([3])), whereas 
the momenta of the incoming and outgoing fermions are collinear with an energy of 0{Mw). 
Also the soft photons couple to collinear fermions with eikonal vertices ~ k^, that do not 
change the polarization of the particle the photon is attached to. Even if these arguments 
would underestimate the effect, the numerical results in the fourth column of table [T] show 
that a 10 per-cent uncertainty on the effect of the cuts on the soft corrections translates to 
an uncertainty of the cross section much below 0.1 fb and is therefore not relevant for the 
desired accuracy needed for the mass measurement. In summary, while a fully differential 
calculation of the N^' ^LO corrections with the EFT method is not currently possible, the mild 
cuts relevant to the W mass measurements together with the smallness of the correction allow 
us to conclude that adding the inclusive N^'^LO cross section estimate to the NLO result is 
sufficient in practice. 

6 Conclusions 

The aim of this paper has been to ascertain that the four-fermion production cross section 
near the P^-pair production threshold can be calculated with an accuracy that matches the 
estimated 6MeV error in the Vl^-mass measurement from an energy scan in e~e^ collisions, 
which requires the consideration of NNLO terms in the Standard Model. The effective field 
theory approach to unstable-particle production makes it possible to investigate the relevant 
part of the Standard Model NNLO diagrams without facing the difficulties of a full calculation, 
which would be overwhelming at present, and to identify a set of parametrically enhanced 
NNLO corrections, all associated with the strong electromagnetic Coulomb attraction of the 
intermediate W bosons. We obtain compact analytic expressions for the total cross section, 
and a result for soft, hard and collinear radiative corrections that extends to the presence of 
any number of Coulomb photon exchanges. The methods used in this calculation may be useful 
in other processes involving unstable particles such as top quarks or new heavy states. 

We find that individual correction terms to the four-fermion cross section are in the 0.3% 
range, leading to shifts of the W mass around 5 MeV. However, the combined effect is somewhat 
smaller. We may thus conclude that the inclusive partonic four-fermion cross section near the 
VF-pair production threshold is known with sufficient precision. We then investigated the 
impact of cuts, showing how to systematically include invariant-mass cuts into the EFT in 
two cases differing by the scaling of the cut with the other parameters of the problem. The 
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cuts applied in the measurement of the M^-pair production cross section at LEP2 loose only 
about 10% of the total cross section. In view of the smallness of the corrections found here, no 
detailed differential cross section calculations beyond NLO appear to be necessary. Our result 
should then be combined with NLO calculations of the four-fermion cross section in the SM in 
the complex mass-scheme [ll^ [T2] as described in appendix [Bj 

The partonic four-fermion cross section must be convoluted with electron structure functions 
that sum large collinear logarithms. As discussed in [13] to achieve the required precision for 
the VF-mass measurement it will be necessary to improve the treatment of initial-state radiation 
that was applied at LEP to account consistently for next-to-leading logarithms. This problem 
is common to all high-precision studies in high-energy e~e~^ collisions. 
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A Renormalisation of the Coulomb potential by hard correc- 
tions 

In this appendix we discuss several technical aspects related to the hard corrections to the 
Coulomb potential used in section [3T6l and their dependence on the renormalisation scheme. The 
calculation of the corrections to the single-Coulomb exchange requires a matching calculation 
where one computes the renormalised W^W~ — > W^W~ NLO scattering amplitude in the 
full theory for (pi -|- p2)^ = 4M^ and compares to the one-loop il^fi^ -^ Q^^^ amplitude 
in the effective theory. Equivalently, the full-theory calculation can be split into contributions 
from different momentum regions. The relevant regions are the hard, potential, soft and semi- 
soft regions. The only contribution that will not be reproduced by diagrams in the EFT is 
that from the hard region, so for the matching calculation it is sufhcient to calculate the hard 
corrections. At leading order in the non-relativistic expansion only the corrections to the single- 
Coulomb exchange diagram contribute. In lA. II we define our renormalisation conventions, the 
hard corrections to the process W^W~ — > W'^W~ are discussed in lA.2[ The relevant results 
in the a{Mz) and G^ input parameter schemes are collected in lA.3l 

A.l Charge renormalisation 

The lowest perturbative scale relevant near the V^-pair production threshold is the W width 
Tw, so we will employ a renormalisation scheme S for the electric charge that is not sensitive 
to smaller scales, in particular not to the light-fermion masses. In practice we will use a{Mz) 
or the Fermi constant G^ as input parameter (see e.g. [31]), but for the moment we will leave 
the renormalisation scheme unspecified. Following the renormalisation conventions of [13] the 
overall one- loop counterterm of the W~^W~ -^ W~^W~ amplitude in the full theory in a given 
charge-renormalisation scheme 5 is given by 

Afounter = [tree] X (2 SZ^ + 26Zw) , (54) 
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since the tree-level single photon exchange diagram (denoted by [tree]) is proportional to e^ = 
(47rQ)^ and has four external VK-legs. 

In the following we write the charge counterterm in a given scheme S as the counterterm 
in the a(0) scheme and a finite scheme-dependent shift 



1 



6Z^ = 5Z°(°) - -Aa^, 
with the explicit expressions for the charge counterterm in the a(0) scheme 



(5Zf(0) 



1 



-SZaa 



S7/1 1 



tAZ 



SZzA 



2 di? ''='=° ^ ^^vf| 



(0) 



Z Ci[j Z 

where the transverse self-energies are defined by the decomposition 



Cii) = ( 5m. - ^ ^Vii') + ^nr(^^) 



(55) 



(56) 



(57) 



For the argument given below we assume that light-fermion masses are used as IR regulators 
in the a(0) scheme, but this will only be used in intermediate steps and the dependence on the 
light-fermion masses drops out in the end. 

A. 2 Hard corrections 

As discussed at the beginning of this appendix we need to calculate the hard corrections to the 
VF"'"VF~ — > W~^W~ amplitude for external momenta directly at threshold. The contributing 
diagrams are of the form of box corrections (figure [T2b) . vertex corrections (figure [T2b) and 
self-energy insertions (figure [T2b.d). The box corrections with a hard loop momentum do not 
contribute at the order we are considering since near threshold they are suppressed by a factor 
V compared to the diagrams with a Coulomb photon. This leaves the vertex and the bubble 
corrections. 





W W 



w w 




Figure 12: Sample diagrams in the full SM contributing to the hard corrections to the WW -^ 
WW subprocess 

Expanding the vertex-correction diagrams of the form of figure 112b in the region where 
the photon and the Ws attached to the vertex sub-loop are potential and the momentum 



Here all conventions are the same as in [34] apart from replacing E — > — D. Note that [34] defines the 
vacuum polarisation TiAAik'^) = ^AA{k^)/k^ which we don't use in the following. 
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running in the vertex loop is hard, one obtains the single-Coulomb exchange diagram with 
an insertion of the one-loop WWj vertex function evaluated with on-shell external momenta. 
The renormalised on-shell vertex function vanishes in the conventional on-shell renormalisation 
scheme where a(0) is used as input parametero Therefore the diagrams with insertion of a one- 
loop unrenormalised WWj vertex are given by the negative of the corresponding counterterm 
in the on-shell renormalisation scheme: 

Avertex = [tree] X {-25^w^) = [tree] x (-2) (dzf^^ + 5Zw + ^SZaa - \—SZza) • (58) 

A similar argument shows that in the hard region the self-energies in figures [T2b/d are 
evaluated at zero external momentum. Since the (unrenormalised) one-loop photon self-energy 
vanishes at zero external momentum the first non-vanishing contributions of diagrams of the 
form 112b comes from expanding the self-energy to first order: 

A, = [tree] x (zH^ '^-'^(fc^))-! = [tree] x {^^^^\^-=}\ , (59) 

where Hj^^"^^ includes all particles except the light fermions. Here it was used that for a hard 
loop momentum the light-fermion masses must be set to zero so the loop integral is scaleless 
and only the heavy particles contribute to the hard part of the self-energy at zero external 
momentum. The unrenormalised photon-Z mixing diagrams give (suppressing the transverse 
projectors) 

A,/,^[tree]x2(.n-(0))^^^ 

.nMeavv(O) 




[tree] x -2 -'—/^ ( "f ) = [t^^^l >< (^^za) ( -f ) , (60) 



since for a potential momentum /c^ <^ Af^. Here the on-shell scheme definition 5Zza = 
— 2n^^(0)/M^ has been used. We also used that massless fermion loops do not contribute to 
the 7 — Z mixing at zero external momenta, as can be seen from the explicit one-loop result 
in [M]. Therefore n^f^^^y(O) = n^^(O). 

Adding the vertex correction, the 7/Z mixing and the counterterm (j54p (where 'tree' is 
again only the photon exchange diagram) one obtains the hard correction in the renormalisation 
scheme S: 

^hard = ^vertex + A^ -^ A^/z + Af^^^te, = [tree] X jf^^j, (61) 

where we defined the correction factor 

6L, = 2{5Zl - 5Zf'^) - 5Zaa + ^^%^U^=o = "Aa^ - ^^^^1.^=0- (62) 

In this expression the derivative of the light-fermion contribution to the photon self-energy 
depends on light-fermion masses used as regulators in the a(0) scheme. This dependence will 



^This relies on 'charge universality' in the standard model, i.e. on the fact that the on-shell electron-photon 
vertex and the on-shell VF-photon vertex receive the same radiative corrections [43] . 
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be cancelled by a similar term in the conversion factor Aa^ for a scheme S that is not sensitive 
to scales below T\y. 

Alternatively, the result (I62p is obtained by considering the individual renormalised contri- 
butions instead of applying the overall counterterm (j54p . In this case 



ALd = A;ertex + A!; + A!;/^. (63) 

Since the on-shell renormalised Z — 7 mixing two-point function vanishes at zero momentum 
we have AJ" /^ = 0. The renormalised one-loop correction to the WW'j vertex vanishes in the 
on-shell scheme so the only non-vanishing contribution comes from the change in the charge- 
counterterm: 

A^ertcx = 2 [tree] x (<5Zf - <^Z-(o)). (64) 

The renormalised self-energy correction to the photon-exchange is given by (I59p and the cor- 
responding counterterm: 

a; = [tree] x I ^^g^ \ .^o " ^Z^^j ■ (65) 

A. 3 Formulas for the a{Mz) and G^ schemes 

We now specialise the result (|62p to the two schemes used in the main text. In the a{Mz) 
scheme the finite shift of the charge counterterm is given by (see e.g. 



^" - Sfc^'^'=° + Ml ■ ^^^^ 

Inserting this definition into (j62p , the sensitivity on the light-fermion masses drops out and one 
obtains the final result for the hard corrections to the Coulomb potential (charge counterterm) 
in the a{Mz) scheme: 

,y, ^ Renfi,(M2) 

In the Gp, scheme the shift in the counterterm is instead given by the correction to muon 
decay, Ar, 

AoP^-Ar- ^Et1(^|, 2'^'- 2^-^^^^°) n^(0)-Ren^(M^) 



bs^ _ lei f^eU^{M^) Ren^ (M| 



s^ 2sl\ M^ Ml 



(69) 



The explicit result for the hard corrections in the ag scheme from ()62p reads 



\ard - -^^ Qj^ \k^=0 

_ ^ni^ely(fe') | , js^ , ,c^ nj^^(o) , ny(o)-Ren^(M^) ^ 
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(71) 



To convert the result from the a^M^) scheme to the G^ scheme we have to add 

Oa{Mz)^Gt, - "hard ~ "hard " J^ ^ 'q^ 1^2=0 

, js^ , ,c^ n^^(o) , n^(o)-Ren^(M^) 

Exphcit expressions for the self-energies appearing in these quantities can be found for instance 
in [Ml. 



B Conversion to the complex-mass scheme 

In the main text we describe the calculation of those N'^' ^LO correction in the EFT framework, 
which are not already included in a NLO calculation in the full Standard Model. The implicit 
assumption is that the NLO diagram corresponding to exchange of a single (Coulomb) photon 
between the W bosons is calculated in the on-shell scheme and defined with a fixed-width 
prescription for the W propagatorsjj Here we discuss the changes that should be made, if this 
diagram is calculated as part of a complete NLO SM calculation in the complex-mass scheme 
defined as in [mill]. 

Since the complex-mass scheme uses the same VK-propagators as the fixed-width scheme, 
the differences arise only from the different renormalisation conventions for a^w and a, and the 
different treatment of H^^'^^ from VF-field renormalisation. The different coupling definitions 
lead to a change in the value of (y 'j^ in the C x H correction (j20p and an additional term 
involving dCf^^acMs ^^ ^^^ NLO Coulomb potential term (f39l) . These compensate exactly for 
the difference between the coupling prefactors a^^a in the single-Coulomb exchange diagrams 
in the complex-mass scheme and the G^ scheme. Thus when the NLO SM calculation is done 
in the complex-mass scheme, we must add 

Aa£i - [Aa£i]cMS (73) 

expanded to order a'^^^a'^ to the N^'^LO calculation performed in the main text. However, at the 
indicated order the correction (j73p actually vanishes. To see this, extract the relevant coupling 
constants from the LO scattering amplitude for four-fermion production A^^' = OewA^ , and 
the amplitude including single-Coulomb exchange, A^^''^'''^'^ = OewoA^'^'^'''^^. The correction 
to the cross section then involves the quantity 



2 Re 



_4(0) f_4(i/2),cA = 2|a„„PRe a^W a(i/2),ci 



(74) 



In the implementation of the complex-mass scheme of [12j the imaginary part of the elec- 
tromagnetic charge arises only through complex masses entering loop integrals in the charge 
counterterm. Therefore the imaginary part of a involves terms of the form a^ Im[log((M^ — 
iMw'rw)/fJ''^)] ~ a^T]y/M\Y ~ a^aew and does not contribute at order Og^a^ in (173]) due to 
the extra factor of T\y/M\y ~ Uew Similarly the difference between a^^ and |ae«)P is of higher 
order, al^{rw/Mwf ~ aL- 



^Since all N^'^LO corrections calculated in the paper are 0{a) corrections to the NLO diagram with a 
single-Coulomb photon, other NLO diagrams are not relevant. 
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A change in the renorniahsation convention that changes the renorniaUsed value of II^^'^) 
affects the calculation of the residue correction, but also the value of c '^^ and the NRQED 
WW'j vertex, which depend on the M^-field renorniahsation convention. The two exactly 
compensate each other as can be seen as follows: a change in the real part of II'^^'^^ changes 
the residue correction by an amount that is directly proportional to Acr£^, as is the term from 

the change of c„'^^ , so the cancellation is obvious. A change in the imaginary part of II'^'^^ 

is more subtle, since it affects the imaginary part of cj, ^^ . We dropped this imaginary part 
in the on-shell scheme, since it does not correspond to cuts involving the four-fermion final 
state, but this is no longer correct if the field renormalisation is complex. In this case one must 
include the imaginary part of c^ ^^ that comes from II*-^'^', in which case the C x H correction 
is no longer proportional to Acr£^. Rather, as in the calculation of the residue correction, one 
obtains the integral (f47l) . We thus see again that the change in H^^'^' is compensated by the 
one in c^ '^^ . A similar argument applies to the change in the WW'j vertex. Hence no further 
correction is needed to convert our result to the case when the full NLO calculation is done in 
the complex-mass scheme. 
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